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Abstract 

In this paper we proposed to use the group of analysis of symmetries of the dynamical system 
to describe the evolution of the Universe. This methods is used in searching for the unknown 
equation of state. It is shown that group of symmetries enforce the form of the equation of state 
for noninteracting scaling multifluids. We showed that symmetries give rise the equation of state 
in the form p = —A + w\p(a) + W2a 13 + and energy density p = A + poia~ 3 ^ 1+w ^ + po2«^ + P03d~ 3 , 
which is commonly used in cosmology. The FRW model filled with scaling fluid (called homological) 
is confronted with the observations of distant type la supernovae. We found the class of model 
parameters admissible by the statistical analysis of SNIa data. We showed that the model with 
scaling fluid fits well to supernovae data. We found that fi mj o — 0.4 and n ~ — 1 (j3 = — 3n), 
which can correspond to (hyper) phantom fluid, and to a high density universe. However if we 
assume prior that f2 mj o = 0.3 then the favoured model is close to concordance ACDM model. Our 
results predict that in the considered model with scaling fluids distant type la supernovae should be 
brighter than in ACDM model, while intermediate distant SNIa should be fainter than in ACDM 
model. We also investigate whether the model with scaling fluid is actually preferred by data over 
ACDM model. As a result we find from the Akaike model selection criterion prefers the model 
with noninteracting scaling fluid. 

PACS numbers: 
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I. INTRODUCTION 



While the structure of space-time is governed by the field equations the physical properties 
of matter, introduced via energy-momentum tensor, is take from other parts of physics. 

Mathematically, one proceeds in an analogous way to that adopted by Newton who 
defined a dual geometric-physical concept of mass-point. In General Relativity one defines 
particles as curves in space-time and ascribes to them physical quantities such as density and 
pressure. The particle of the rest-mass m is defined as a future directed curve 7 : / — > M in 
a space-time M, such that #(7*, 7*) = —m 2 . The tangent vector 7* is the energy-momentum 
vector of the particle. Then, a particle flow of rest-mass m is defined to be a pair (P, 77) 
where 77 is a function rj : M — > [0, 00) called the world density, and P is the energy-momentum 
vector field P: M —>■ T(M) such that each integral curve of P is a particle of rest-mass m. 
The energy-momentum tensor of particle flow (P, 77) is tensor T = 77P ® P; and pressure p, 
understood as a function of M, enters it through P (for detailed formalism see y], Mc 
Crea [3] noticed many years ago that since matter is taken from outside of general relativity, 
there are no a priori reasons why p should be non-negative. The only way how p should 
be interpreted is through effects which itproduces in a model. For example, p < may 
be interpreted as a quintessence matter 4J or dark energy due to the present Universe is 
accelerating P, |(| . 

Lie groups of symmetry play an important role in searching for new solutions but their 
original field of fruitful applications remained hidden in the literature. This problem was 
discussed by Stephani j^] in his introduction to the basic monograph on the application of 
symmetry groups to general relativity. But this method is still obscure because many people 



who used them implicitly simply were not aware of their existence 

The equation of state also plays an important role in general relativity. As it is well known, 
Einstein's field equations, together with the Bianchi identities, form an undetermined system 
of equations and one more equation, equation of state must be added. Einstein's equations 
completed by the equation of state satisfy certain symmetries which has structure of a group, 



called a symmetry group (see It is interesting that this procedure can be reversed 

. Collins applied the method of deriving form of the equation of state from symmetry 
to the case of classical and relativistic stars to obtain the physically realistic equations of 
141 ] . This illustrate-to use Collin's expressions-a subliminal role of mathematics 
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in physical considerations which consists in suggesting correct physics on the ground of a 
logical beauty (such as, for example, symmetry principles). 

There are several interesting papers where the authors study symmetry transformations 
under which the Einstein equation are invariant. They relate symmetry and inflation ^, 
0, Q] however the Lie group of symmetries is not used in analysis of the FRW equation. 
The investigation of Lie symmetries of General Relativity and cosmology and the origin of 
dynamical systems in the context of casual viscous fluid and the cosmology with variable 
constants can also be found interesting. For a cosmological models with bulk viscosity idea 
of renormalisation group is also presented to study some scaling properties of the model 
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Because the Einstein equation constitute very complicated system of nonlinear partial 
differential equation it is not easy to obtain an exact solution. However, one can assume 
self- similarity symmetry and then due to this simplified assumption Einstein field equation 
with spherical symmetry can be reduced to the form ordinary differential equations i.e. a 
dynamical system. Therefore the assumption of self-similarity symetry is very useful in the 
problem of finding exact solutions in nonhomogeneous case (see [3, H])- The self-similar 
solutions play especially important roles in the cosmological applications or description of 
gravitational collapse. It is the consequence of the fact that they describe asymptotic be- 
haviour of more general solutions. The idea that spherical symmetric fluctuations might 
naturally evolve toward self-similar form is kaowa as as the similarity hypothesis Q. 

A very interesting result was founded by Harada and Maeda [121] . They proved numer- 
ically that self-similar solution acts as an attractor in the phase space for spherically sym- 
metric collapse of perfect fluid. In the context of gravitational collapse self-similar solution 
has been also studied from the points of view of critical phenomena j^j]. 

The self-similar solutions also play an important role in the contemporary context of dark 
energy. The scaling solutions in the context of "cosmic coincidence conundrum" were also 
considered in j^j]. It is demonstrated in many papers that there exists a stationary sca ling 
solution in this problem which is representing a stable attractor at late times 13, |2jJ, |22j, |23j, 
I24 ] . Among the class of solutions, the global atractor solutions are of the greatest interest 
from the physical point of view. They usually constitute scaling type of solutions which can 
explain why energy densities of dark energy and matter at present epoch are of the same 



order of magnitude 
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Our main idea is to use the concept of self-similarity to the derivation of the form of 
equation of state. They reproduce themselves as the scale changes. This method olog y was 
already used by Ellis and Buchert in context of problem of gravitational entropy 29( . 

To determine the structure and the dynamics of astrophysical systems as well as the 
Universe, the equation of state is usually necessary. Let us consider, for example, the 
structure of a spherical neutron star. Then if the pressure p is known as a function of the 
density p, we can determine the gravitational mass M and the radius R of the star as a 
function of the central density p c by solving the Oppenheimer-Volkoff equation. This means 
that we can determine in principle the mass-radius relation M(R) theoretically. In similar 
way if the dependence between total pressure and total energy density of the Universe is 
postulated then we can integrate the Friedmann equation |3Cj. However, the equation of 
state relevant to the neutron star or the Universe is not established yet, although it may 
be determined in the near future. Let us note that for the Universe the form of equation 
of state is not theoretically known at present but from observations of SNIa data it may be 
reconstructed in the SNAP3. Q]. 

In our paper the form of the equation of state is not postulated a priori but it is derived 
from the existence of scaling solutions. The existence of scaling solutions is equivalent the 
existence of a corresponding symmetry group (similarity group) of the basic differential 
equations. It is just discussed by us the type of symmetry of dynamical equations. Our 
philosophy is to assume the existence of scaling solutions instead of postulated a priori any 
prescribed form of the equation of state. 

We use Lie symmetries to determine the dependence of energy density on the scale factor. 
Of course if we assume the exact form of the equation of state for example for phantoms 
then it is easy to obtain the exact form of p e ff{a) for the adiabatic condition. 

In physical applications the important role is played by the scaling solutions, they are also 
called self-similar, homological etc. The physical properties of such solutions are important 
in many branches of physics 1£, 0, 2^]. They can be found in the exact form or can be 
analysed by using the dynamical system methods. Then invariants of self-similarity groups 
are a good choice for a variable to parameterise a system. 

A natural question is whether the cosmological model filled by scaling fluid is actually 
preferred by the data over concordance ACDM model. This problem can be investigated 
by using the Akaike (AIC) and Bayesian (BIC) informative criterion. We find that AIC 
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favoured FRW model with scaling fluid relative to ACDM model. 



In the present paper we pursue Collins' way of thinking |13l . 114 ] applying it to the field of 
cosmology. In section 2 we outlined the theory of symmetry groups of differential equations, 
then we showed that cosmological equations, together with a suitable equation of state 
admit a certain Lie group of symmetries or, vice versa, the invariance of equations with 
respect to a given symmetry group singles out corresponding equation of state. In sections 
3-4 we showed this for Friedmann models with matter in the form of perfect fluid and 
effective cosmological constant, respectively. In section 5 and 6 we confronted data obtained 
theoretically from symmetry equation of state with recently available SNIa data. In section 
7 we briefly commented on the obtained results. 

II. SYMMETRY GROUP OF A SYSTEM OF DIFFERENTIAL EQUATIONS. 



,u m ), n + 



Let us consider, in a Euclidean space E (x, u), x = (x 1 , . . . , x n ), u = (u 1 , 
m = N, the system of differential equations 

du a 

p a = f a (x,u), « = (l,...,m), P a = ^ (!) 
and point-point transformations 



. x = x(x, u), t 
T : I (2) 



u a = u a (x, u), 

which map each solution of system into a solution of the same system. T is a Lie group. 
(All these considerations may be easily generalised to the case when differential equations 
are defined on a A^-dimensional differential manifold; see [3, 

Let M: u = u a (x) be a solution of (P). M defines a submanifold in E N ; if M = T(M) 
we have 

— \ x = x(x, u(x)), / 
M:\ , , ' (3) 

I u a = u a (x, u(x)), 

The derivatives p a = and p a = satisfy the condition 

p a Dx = Du a , (4) 

where D = J| +P a -^- By solving (jlj) we obtain p a = p a (x, u,p). By joining these solutions 
to transformations T, one gets a new set of transformations T which is called "extension of 
T to the first derivatives" . 
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If the infinitesimal operator of the group T has the form 



X = ^ u )^ + V a (x,u)^ (5) 



then that the group T is 



X = X+(Drf- V PD£)-^ (6) 

Function T defined on E (x,u,p), which are preserved under transformations T, are called 
invariants of the group T; they satisfy differential equations 

X{F) = (7) 

and, vice versa, functions satisfying (J7J) are invariants of the group generated by X. 
Now, we may apply equations (J2J) to our original equations ((TJ) to obtain 

X(p a -r(x,u)) = (8) 

Putting (JOJ) into (jSJ) and taking into account p a = f a (x, u), we get 

drf + drf _ d Afa _ = ^d_r + dJ1 ^ ^ 

dx dui 3 dx du@ dx du 13 ' 



III. THE GROUP SYMMETRY OF EINSTEIN EQUATIONS WITH THE PER- 
FECT FLUID ENERGY-MOMENTUM TENSOR. 

The problem of group properties of vacuum Einstein equations has been investigated 
by Ibragimov (1983) 33j. In this case we obtain one parameter group of rescalings of 
the metric gij — > agij together with infinite group of coordinate transformations from the 
maximal group of point transformations admitted by vacuum Einstein equations Rik = 
0; i, k — 0, 1, 2, 3. We shall investigate symmetries of Einstein equations: 

R ik = T ik -^Tg ik i,k = 0,1,2,3 (10) 

We assume that Tik is the perfect fluid energy-momentum tensor. We choose the comoving 
coordinate system i.e. U{ = 5 0i . 
Then we have: 

T ik = (p + p)Soi6 ok - pgik 

T = p-3p, (11) 
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where p = p(p). 

We assume that Einstein equations (fTUJl with energy-momentum tensor (fTTj) possesses the 
same symmetries as corresponding vacuum equations. We shall look the form of equation of 
state enforced by symmetries of the vacuum Einstein's equations. Infinitesimal operator of 
the Lie group symmetry acting in the set of Einstein equations can be written in the form: 

* = rtx)A + *(»,p)A + S(l ,,, rt ^. (12) 

where 9 can be derived from the condition for operator (j!2j) to be admissible operator of 
equation (JOJ). This condition is the following: 

X 2 [R ik - T ik + -Tg ik ] = 0. (13) 

The operator X 2 is the extension of X to the second derivatives: 

X 2 = X 2 + 9— 
op 

X 2 = X + iijk^ 1" 

°9ij,k 09ij,kl 

X 2 [R lk ] = (14) 

It was demonstrated that symmetries of the vacuum Einstein equations enforce the equation 
of state p = const + wp (w = const) if we postulate that the source of gravity is given in the 
form of perfect fluid |3Ja] . 

IV. SYMMETRY GROUP OF FRIEDMANN EQUATIONS 

Equation © gives us the conditions of the existence of the operator symmetry X for the 
system (JIJ. Algebra of these operators characterises the symmetry of the group T symmetry 
of system (JTJ). 

In the following part we shall consider symmetry of autonomous dynamical systems for 
which x a = f a (x) (a dot denotes the differentiation with respect to the parameter u), 
^7 = and infinitesimal transformations are generated by X = £(x) + rj a [u)-^. In this 
case (JHJ) simplify to the form 

du 13 ^ dx^ ^ du 13 ^ ^ 
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The Friedmann equations may be written in the form of autonomous dynamical system on 
the plane (a, p) 



da / pa 2 Aa 2 , 

a = Tt = y—- k+ — = f{a > p) 



p= d {t = -- a f-T- k + ^t(p+pM) = AM (16) 

where a- is a scale factor, p- energy density, k- curvature index, A- cosmological constant. 

The set of basic dynamical equations is called the Friedmann dynamical system. If we 
assume that the operator of infinitesimal transformation of symmetry for system (JTfij) is of 
the form 

x =^m + ^Ta + ^r P < 17) 

then equations (fT5|) can be reduced to 

Let us concentrate now (without lose of degree of generalisation) on the special case of the 
FRW, namely the flat model without cosmological constant. Then equations (fT7)J) reduces 
to 

^ = aE 
dt 

d -§ = -\^ + P ^E)) (19) 

where E = \/~p/3 is redefined energy density, p is a pressure as an unknown function of scale 
factor. Let us assume that equations (|T9*j) admit the special symmetry of the operator for 
quasi- homology, i.e. 

x =«4 + " I(a) l + " 2(£) 4 (2o) 

Then from equations (fTH|) we obtain 

d^_<m = ^ + rj\E)_ 
da dt a E 



drf(E) d£ _ 1 ( m \ i „\ d ran , d P 

dE dt ~ 3E 2 +p(a,E) 



{^Ya + {QE + ^ )r]2{E) ) (22) 



S 



^From (|21|) we have that 

r] 1 = Aa, 7] 2 = DE, f (t) = Bt (23) 

where A, D, B are constants which satisfying the constraint condition: 

D = -B (24) 

After substitution (J2*3~j) and (J2H) into the equation (J2*2~j) we obtain partial differential equation 
determining p = p(a, E) 

2D P =Aa d £ + DE^ (25) 

The solutions of equation (f2*K|) determine the form of the equation of state required by the 
condition of quasi-similarity (quasi- homology). The condition for the existence of quasi- 
homological type of symmetry reduces to corresponding condition for existence of homolog- 
ical symmetry and then p must satisfy (|25|). 

Let us consider now some special solutions of (|25|). If p = p(E) only then p oc E 2 = p, 
while if p = p(a) only then p(a) oc . It is easy to find the general solution of (}25|) by 
using the standard characteristic methods. Then we obtain 



E a, a 



ly'— J =0 ( 26 ) 

is the solution of eq. (|25J). where G C 1 , and Ji, J2 are its invariants. 

There is equivalent form of the solution (j2*U|) in which the pressure can be given in the 
exact form, namely 

or g{x) oc x where g is any function of class C 1 . If we substitute g(x) = const in eq. (|27|). 
then mentioned before special cases of p can be simply recovered. Therefore eq. (J27j) gives 
the most general form of the equation of state for the Universe to be self invariant or 



homological. In this case we can find a strict analogy to stars |l4j. One should notice 



that the formulae (}2Tj) contains for example the following form of the equation of state 

2D 

p = -A + Cip(a) + C 2 a^ (28) 
which is crucial for our further investigations. 
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We can find that eq. (|28|) describes the pressure of noninteracting multifluids components. 
Then from the conservation condition postulating for each component we obtain p = A + 
Poia~ 3( - 1+Cl ^ + po3 a ~ 3 + Po20^, where A, p 01 , p 02 , P03 are constants. Note that in the 
special case of W = — 3 equation (J28)l and conservation condition give rise to new type of 
contribution — (3pm In a)a~ 3 instead po2<2~. It corresponds to MOND phase squeezing in 
FRW scenario 36[. To obtain the form of equation of state (|28|) it is sufficiently to choose 
g(x) = —3 + C% + C2X. Then from equations (|2*3*j) the finite transformations can be simply 
obtain for this case, namely 

— At 

a — > a = ae 
E -> E = Ee Dr 

t^t = te~ D \ rel 1 (29) 

2D 

One can check that the invariant a ^"p is really preserved under the action of the homolog- 
ical transformation, i.e. 

a a p(a,E) = a a p(a,E) (30) 

The equation of state given by ()27|) assumes very general form. Their existence can be 
closely associated with the real equation of state for quintessence matter. Let us note that 
equation ([2~S|) contains not only standard term of the equation of state p oc p, but also 
the cosmological constant with a constant part and a part depending on the scale factor a. 
In the special case if D = —A we recover the form of the effective cosmological constant 
distinguished by a maximum of the amplitude for the tunnelling Gamov process j^]. In 
the general case in the parameterisation of the cosmological term it is usually assumed the 
oower law form of A(a) dependence which is also suggested from symmetry consideration 

Note that homological transformations are strictly related to dimensional analysis. The 
existence of such kind symmetry means that equation can be simply expressed in the di- 
mensionless form. Therefore for each component of the fluid one can define dimensionless 
density parameters which contain basic Friedmann equation. 
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V. EQUATION OF STATE PARAMETER FROM DISTANT SUPERNOVAE OB- 
SERVATIONS. 



In this section we confront FRW cosmological model filled with homological fluid (which 
equation of state was previously suggested from symmetry arguments) with observations of 
distant supernovae type la. These interpreted observations in the framework FRW model 
indicate that our universe is presently accelerating due to unknown form of matter called 



dark energy with negative pressure |39|, |40 
which was tested from SNIa observations 
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.There are different candidates for dark energy 

"nn 



42,143 



4j,|44|4£| 



The most popular candidate is cosmological constant which can be treated as a per- 
fect fluid for which p = —p and p = A. Although cosmological constant explains ob- 
servations of SNIa we don't understand why the observed value of cosmological constant 



%-kG 



* 10 GeV ) is so small in comparison with its natural theoretical expectations from 
quantum field theory ~ 10 71 GeV^ 4 . According to this theorythe energy-momentum tensor 
of the vacuum is nonvanishing and < >=< p > g^ iv . Therefore the observed cosmologi- 
cal constant term is A = A + 8ttG < p > where A is the "bare" cosmological constant. Our 
naive intuitions is that < p >~ Mpj while upper bound on the present value of A (referred 
as A ) can be given in term of Hubble's constant H 

|Ao|/8ttG < 10- 29 g/cm 3 ~ l(T 47 Gel/ 4 

The idea that cosmological constant is the sum of two terms A and 8irGp vac is proposed as 
a method of explaining the observed small value of A which might had been large in the 
early universe. The variable A models (or a decaying vacuum energy density) are described 
in terms of two fluids mixture just like it was predicted from similarity analysis. In the 
majority of the paper's calculations py depends only on the cosmological time through the 



scale factor py ~ a 2 . The expression p 
made in the spirit of quantum gravity 



a) is obtained by using dimensional arguments 
| . We assume that we have a decaying vacuum 
medium parameterised by scale factor in the power of law plus a noninteracting cosmological 
constant term A. Additionally, we usually have baryonic matter satisfying the equation of 



state p = 0. Hence the FRW equation in general case can be rearranged to the form giving 



11 



the Hubble function H(z) = a/ a 



H$ \a J \a J \a J \a ) 

For a = a (the present value of scale factor) we obtain the following constraint 



^m,0 + ^C,0 + ^fc,0 + &w,0 + ^A,0 — 1- (32) 

Because of the constraint H 2 > the space admissible for the motion should be restricted 
to the region at which p e s(a) > 0. If some negative component of effective energy density 
appears in p e &{a) then it can never dominate at the early universe and at the late times. 
If negative energy density scaling like a~ m appears in such a way that 3H 2 = Q m fia~ m — 
fLoa~ n + Qa,o (where n > m) then we have the case of bouncing cosmological models 



48t |49(. In this model the relation p e g(a) > is of course satisfied. In general the 
violation of the strong energy condition (SEC) is a necessary (but not sufficient) condition for 
the bounce to appear. The bouncing models can be characterised by the minimal condition 
under which the present Universe arises from a bounce with a prev ious collapse phase (note 
the reflectional symmetry of the basic equation H —> —H) j3J3- I n or der if we have such 
a situation that n < m, then negative term — f2 nj oa _n cannot dominate late evolution and 
therefore the scale factor should be bounded and an inverted bounce appears. 

If the energy density is very large then quantum gravity corrections are important at both 
the big-bang and big-rip singularities. The account of quantum effects leads to avoid not 
o„ ly the initial S i„ gular it y but a, S o escape from the futute 8 in gU .a r it y QHQ. 

The idea of bounce in FRW cosmologies appeared in Tolman's monograph devoted to 
cosmology 0]. This idea was strictly connected with oscillating models 5fil. I57I Issl] . At 




present oscillating models play an important role in the brane cosmology [59|, |6Cj • The FRW 
universe undergoing a bounce instead of the big-bang is also an appealing idea in the context 
of quantum cosmology j^lj . The attractiveness of bouncing models comes from the fact that 
they have no horizon problem and they explain quantum ori gin of structures in the Universe 



47L l48| characterised the bouncing 



Molina-Paris and Visser and later Tippett 
models by the minimal condition under which the present universe arises from a bounce from 
the previous collapse phase (the Tolman wormhole is different name for denoting such a type 
of evolution). The violation of strong energy condition (SEC) is in general a necessary (but 
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not sufficient) condition for bounce to appear. For closed models it is sufficient condition 
and none of other energy condition need to be violated (like null energy condition (NEC): 
p + p > 0, week energy condition (WEC): p > and p + p > 0, dominant energy condition 
(DEC): p > and p ±p > energy conditions can be satisfied). 

We can find necessary and sufficient conditions for an evolutional path with a bounce by 
analysing dynamics on the phase plane (a, a), where a is the scale factor and dot denotes 
differentiation with respect to cosmological time. We understand the bounce as in 0,2] 5 
namely there must be some moment say t = bounce i n evolution of the universe at which the 
size of the universe has a minimum, dbounce = and a > 0. This weak inequality a > is 
enough for giving domains in the phase space occupied by trajectories with the bounce. 

We understand the inverted bounce as the situation in which in some moment t = bounce 
in evolution of the universe the size of the universe has a maximum, dbounce = and dbounce < 
0. 

For our aim we considered some version above model in which for simplicity we put Q w $. 
If we rewrite this equation in z variable we obtain for the flat model that 



H(z) = H ^n mfl (l + zf + ft Aj0 + (1 - fi m , - fiA,o)(l + z)\ (33) 

where 7 = 3n. Our further task is be to confront the "homological gas" or variable A fluid 
with SNIa data and for this purpose we calculate the luminosity distance in a standard way 

r dz 

d ^ = ^hm (34) 

Further in this section we will use the flat FRW model since the evidence for this case is 
very strong in the light of WMAP data Therefore while talking about model testing 
we actually mean the estimation of both Q\,o and 7 parameters for the best fitted flat model 
with homological gas. Specifically we have tested this cosmology with prior assumption 
^m,o — 0.3 and fi m = 0.05 as well as a model without prior assumption on r2 m> o- 
To proceeded with fitting SNIa data we need magnitude-redshift relation 

m(z,M, fi m , ; 7) = M + 5 log 10 D L (z, fi m , ; ^a,o, 7) (35) 

where: 

D L (z, n mfi ; tt Xfi , 7) = H d L (z, H , Q rnfi ; fl A , 7) (36) 
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is the luminosity distance with H$ factored out, so that marginalisation over the intercept 



M — M — h log 10 H + 25 (37) 

leads actually to joint marginalisation Hq and M (M being the absolute magnitude of SNIa). 
Then we can obtain the best fit model minimising the function \ 2 

I 'theor ^,obs\2 

x 2 = E ( ' ~2 ' (») 



where their sum is over the SNIa sample and <7j denote the (full) statistical error of magnitude 
determination. This is illustrated by figures of residuals (with respect to Einstein-de Sitter 
model). One of the advantages of residual plots is that the intercept of the m-z curve gets 
cancelled. The assumption that the intercept is the same for different cosmological models is 
legitimate since Ai is actually determined from the low-redshift part of the Hubble diagram 
which should be linear in all realistic cosmologies. 

However, the best-fit values alone are not relevant if not supplemented with the confi- 
dence levels for the parameters. Therefore, we performed the estimation of model parameters 
using the minimisation procedure, based on the likelihood function. We assume that super- 
novae measurements came with uncorrelated Gaussian errors and in this case the likelihood 
function C could be determined from chi-square statistic £ oc exp(— x 2 /2) while probability 
density function of cosmological parameters is derived from Bayes's theorem Q|. 

Therefore, we supplement our analysis with confidence intervals in the (^a,0;7) plane by 
calculating the marginal probability density functions 

7>(fi A ,o, 7) oc J exp(- x 2 (fi m ,o, Oa,o, l,M)/2) dM (39) 

with f2 m fixed (fi m ,o = 0.05, 0.3) and 

P(f] A ,o,7)oc [ [ exp(- X 2 (tt mfi ,tt Xfi , 1 ,M)/2)dtt mfi dM (40) 



without fixed f2 m ,0) respectively (a proportionality sign equals up to the normalisation con- 
stant). In order to complete the picture we have also derived one-dimensional probability 
distribution functions for Q obtained from joint marginalisation over 7 and ^a,o- The max- 
imum value of such a PDF informs us about the most probable value of Q (supported by 
supernovae data) within the full class of homological models. 
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VI. SAMPLES USED 



Supernovae surveys (published data) have already five years long history. Beginning with 
first published samples p, |y] other data sets have been produced either by correcting original 
samples for systematics or by supplementing them with new supernovae (or both). 

Because original Perlmutter et al. and Riess et al. samples j^, 0] were completed seven 
years ago, presently the newer supernovae observations are used. Knop et al. 3| have 
reexamined the Perlmutter sample with the host-galaxy extinction correctly applied. They 
chose from the Perlmutter sample these supernovae which were the more securely spectrally 
identified as type la and have reasonable colour measurements. They also included eleven 
new high redshift supernovae and a well known sample with low redshift supernovae. 

We have decided to test our model using this new sample of supernovae. They make 
possible to distinguish a few subsets of supernovae from this sample. We consider two of 
them. The first is a subset of 58 supernovae with corrected extinction (Knop subsample 6; 
hereafter K6) and the second is that of 54 low extinction supernovae (Knop subsample 3; 
hereafter K3). Samples C and K3 are similarly constructed as containing only low extinction 
supernovae. The advantage of the Knop sample is that Knop's discussion of extinction 
correction was very careful and as a result his sample has extinction correctly applied. 

Another sample was presented by Tonry et al. [3] who collected a large number of 
supernovae data published by different authors and added eight new high redshift SN la. 
This sample of 230 SN la was re-calibrated with a consistent zero point. Wherever possible 
the extinction estimates and distance fitting were recalculated. Unfortunately, one was not 



able to do so for the full sample 
further improved by Barris et al. 



This sample was 



for details see Table 8 in Ref. 

66| who added 23 high redshift supernovae including 15 
at z > 0.7 thus doubling the published record of objects at these redshifts. Tonry et al. and 
Barris et al. presented the data of redshifts and luminosity distances for their supernovae 
sample. Therefore, Eqs. (|33j) and (|3*7jl should be modified appropriately [z2| 

m - M = 5 log 10 ( J D L ) T onr y - 5 log 10 65 + 25 (41) 

and 

M = -51og 10 H + 25. (42) 
For the Hubble constant Hq = 65 km s -1 Mpc -1 one gets Ai = 15.935. 
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Recently Riess et al. 0| significantly improved their former group sample. They discov- 
ered 16 new type la Supernovae. It should be noted that 6 of these objects have z > 1.25 
(out of total number of 7 object with so high redshifts). Moreover, they compiled a set of 
previously observed SNIa relying on large, published samples, whenever possible, to reduce 
systematic errors from differences in calibrations. The full Riess sample contains 186 SNIa 
("Silver" sample). On the base of quality of the spectroscopic and photometric record for 
individual Supernovae, they also selected more restricted "Gold" sample of 157 Supernovae. 

Riess et al.'s sample has been used by many researchers as a standard dataset. However, 
for the sake of comparison and illustration we analysed also earlier Knop j 



sample of 

supernovae. This seems to be useful because, as pointed out in the literature, studies 



performed on different SNIa samples often gave different results (see for example (42 , 



44. 



VII. CONSTRAINING EQUATION OF STATE FROM DISTANT SUPERNOVAE 

In order to test the model we calculate the best fit with minimum y 2 as well as we esti- 
ma te the model p_ Usi „ g the m a* imUm hhehhood method Q. Fo r hoth statistical 
methods we took the parameter n in the interval [—3.33, 2], f2 mj o in the interval [0, 1], Qc,o 
in the interval [—1, 1], while an interval for is obtained from the equation (|32j). 

We have tested the models in three different classes of models. At first we analysed the 
data without any prior assumption about fi m ,o- I n the second class we assumed that (2) 
Q m o = 0.05 while in the last class (3) we assumed that f2 m = 0.3. 

The second class was chosen as a representative of the standard knowledge of fi m (bary- 



onic plus dark matter in galactic halos Q] . In the last class we have incorporated (at the level 
of f2 nij o) the prior knowledge about baryonic content of the Universe (as inferred from the 
BBN considerations). Hence this class is representative of the models in which non matter 
component is responsible both for dark matter in halos and its diffuse part (dark energy). 

It is interesting to compare this results with these obtained for the model with A = 0. 
Such model is equivalent to Cardassian model 7^ . Because this mo- 
alternative to ACDM, it was immediately verified by SNIa observations [42], 
now it is possible to perform more precise test of the models using the new Riess sample. 

At first we have decided to test the our model (both for vanishing and not vanishing 
cosmological constant) using Knop et al. 0] and Riess et al. |3| samples of supernovae. 



was proposed as an 



74 



However 
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In all samples we marginalise over parameter Ai. It means that the Hubble parameter is 
fitted from observations, too. 

The results of two fitting procedures performed on different samples and with different 
prior assumptions concerning the cosmological models are presented in (Table . The de- 
tailed results of our analysis for the flat model with vanishing cosmological constant are 
summarised in Table ITT1 This tables refers both to the \ 2 (best fit) and results from maxi- 
mum likelihood method. In both cases we obtained different values of M. for each analysed 
sample. It should be noted that values obtained in both methods are different but differences 
are much smaller for Riess et al. sample than in the case of older Knop et al. sample. 

At first we analysed full model with not vanishing cosmological constant. For example, 
with in the with the maximum likelihood method, first class of models gives values of the 
parameter: Q m fi = 0.38, flc,o — 0.64, Q\ t o = 0.72, n = —2.10, for sample K6, while 
Q mfi = 0.27, Q c ,o = 0.01, Q A)0 = 0.75, n = 0.00, for sample K3. 

With the new Riess sample we obtain: f2 mj0 = 0.46, Qc,o = 0.40, Q\ : q = 0.12, n = —1.10, 
for Silver sample , while f2 mj o = 0.44, Qc,o = 0.34, f2A,o = 0.20, n = —0.90, for Gold Sample. 

This result mean that for Knop sample flc,o part play only marginal role and the model is 
very close to ACDM model. However with the new Riess sample we obtain that fi m — 0.4 
with n ~ — 1 which correspond to (hyper) phantom model. However if we assume Q m fi = 0.3 
from independent extragalactic estimation jj], then £lc,o is small and n is close to zero, so 
ACDM model is favoured. 

On Fig. 1 we present residuals plots with respect to the Einstein-de Sitter model for 
ACDM and our model. We observed that distant SNIa should be brighter (in our model) 
than in the ACDM model. What is interesting is that the Hubble diagram for our model 
(both for vanishing and non vanishing A) intersects the corresponding ACDM diagram. In 
such a way, the supernovae on intermediate distant are fainter then expected in the ACDM 
model. Please note that the Riess et al. sample contain only very few supernovae in the 
intermediate distance, so this prediction should be tested with future supernovae sample 
when more type la supernovae measurement in the intermediate distances will be available. 

It should be noted that knowing the best-fit values alone have not enough scientific 
relevance, if confidence levels for parameter intervals are not presented, too. Using the 
minimisation procedure, based on the likelihood method we also carry out the errors of the 
model parameters estimation. On the confidence level 68.3% we present parameter values for 
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samples K6 and K3, Silver and Gold (Table ITTTj) . The density distribution (one dimensional 
PDF) for model parameters obtained by marginalisation over remaining parameters of the 
model are presented in Fig. 2-5. Additionally the on the confidence level 68.3% and 95.4% 
is marked on the figures. 

Please note that both positive and negative values of Qc,o are formally possible. It is the 
reason why we received on the Figs 4 and 5 bimodal distributions. Figs 6 and 7 explain this 
situation in more details. These figures present the confidence levels on the plane (f2 m ,o, ^c,o) 
(Fig. 6) and (Qc,o, n ) (Fig- 7) minimised over remaining model parameter. 

Fig. 7 is more complicated than in the case of the model with Q\ = (see Fig. 4 in 



Rcf. 



42j where the maximum likelihood procedure suggests that n should be negative and 



consequently fi m is greater than 0.3). Now we obtain also the possibility that n > and 
^m,o > 0.3 as a result of presence both Qc,o an d terms. 

The density distribution for model parameters for the model with fixed fl m ,o — 0.3 is 
presented on Figs 8-10. Those figures confirmed that in this case the model is close to 
ACDM model. 

For the model with vanishing cosmological constant Qa = the errors of the model 
parameters estimation is presented on Table IIVI For this model, using the maximum likeli- 
hood method, we obtained with the sample K6 that il m ,o = 0.52tQ Qg and n = — 3.33 +2,00 on 
the confidence level 68.3%. In turn, for sample K3 we obtained that fi m0 = 0.481q53 and 
n = — 0.40^24 on the confidence level 68.3%. With the new Riess et al. sample we obtain 
^m,o = 0.511q!o5 and n = — 1.60^1$ with the Silver sample, while with the Gold sample 
^m,o = 0.51^05 and n = — 1. 231^23 ■ The best fit procedure also suggests that n should be 
negative and consequently il m is greater than 0.3. 

One can see that that result obtained with maximum likelihood method, for f2 mj o and 
flc,o are similar for all samples, however with new Riess et al.'s sample errors in parameter 
estimation significantly decreased. Please note that with assumption that A = the model 



is equivalent to the Cardassian models. We confirmed our previous results 42] on the new 
Riess et al.'s sample. The observations favour the high density universe with n negative. 
On the other hand, if we assume fi m> o = 0.3 then n ~ which correspond to ACDM model, 
like for scaling multifluids model with non vanishing A. 

For the model with f2 mj0 = 0.3 we obtained for sample K6 n = —0.13 with <y{n) = 0.23. In 
turn, for sample K3 we obtained n = —0.20, cr(n) = 0.17. With the new Riess et al. sample 
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we obtained n = —0.07 with a(n) = 0.10 for Silver sample, while n = —0.03, cr(n) = 0.10 
for Gold sample. On can see that result obtained with both samples are similar but with 
the Riess et al. sample errors in estimation of the parameters significantly decreased. These 
results mean that if we assume fi m = 0.3 then n ~ which correspond to ACDM model. 

In this way, it is crucial to determine which combination of parameters give the preferred 
fit to data. This is the statistical problem of model selection The problem is usually 
the elimination of parameters which play insufficient role in improving the fit data available. 
Important role in this area plays especially the Akaike information criterion AIC This 
criterion is defined as 

AIC = -2\nC + 2k (43) 

where £ is the maximum likelihood and k is the number of the parameter of the model. The 
best model is the model which minimises the AIC. 

The Bayesian information criterion BIC introduced by Schwarz j3] is defined as: 

BIC = -2\nC + d\nN (44) 

where N is the number of data points used in the fit. While AIC tends to favour models with 
large number of parameters, the BIC the more strongly penalises them, so BIC provides the 
useful approximation to full evidence in the case of no prior on the set of model parameters 

Q- 

The effectiveness of using these criteria in the current cosmological applications has been 
recently demonstrated by Liddle [3] Please note that both information criteria values have 
no absolute sense and only the relative values between different models are physically inter- 
esting. For the BIC a difference of 2 is treated as a positive evidence (6 as a strong evidence) 
against the model with larger value of BIC [sol . Isi| . 

The AIC for the models under consideration is presented in the (Table |^) while for BIC 
in the (Table IVTj) . The BIC information criterion do not show significant differences between 
scaling multifluid model with vanishing A and ACDM model. The model which minimising 
AIC is scaling multifluids model with the vanishing cosmological constant (equivalent to 
the Cardassian model). It means that, from the statistical point of view, if we compare 
scaling multifluid models with vanishing and non-vanishing A, the extra term with A does 
not improve significantly the quality of the fit. 
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Basing on AIC information criterion, it is clear that the scaling multifluid model with 
vanishing A fits better to data then the ACDM model. However scaling multifluid model 
with vanishing A indicate high density universe, with Q m )0 — 0.5 which seems to be too high 
with comparison with present extragalactic data Scaling multifluid model with non- 
vanishing A allow density of the universe Q m ,o close to 0.3 (Fig. 7) while ACDM model predict 
O m0 — 0.3. It clearly shows that more precise measurements of £l m ,o from independent 
observations is necessary to final discriminate between the ACDM and scaling multifluid 
models. 



VIII. CONCLUSION 



The Supernova Cosmology Project and the High-Z-Supernova Search reported of their 
observations of type la supernovae and suggest that the expansion of the Universe is still 
accelerating due to the presence of unknown form of matter called dark energy. For the 
accelerating Universe the equation of state parameter w = - for dark must satisfy w < — g. 
The cosmological constant A is arguable, but at the some time the simplest candidate for 
this dark energy, although it is well known, that predictions for its value are many orders 
of magnitude off from the observationally acceptable value. The introducing of A = A (a) or 
quintessence is usually proposed as a possible solution to avoid the cosmological constant 
problem. 

In presented paper the idea of derivation the form of equation of state from symmetries 
of self-similarity of the FRW dynamics is considered. We have shown that these symmetries 
enforce appropriate equation of state used commonly in cosmology. The property of ho- 
mology is very important in astrophysics, when main sequence in the Hertzsprung-Russell 
diagram can be reconstructed from the invariants of homology transformations. Moreover 
from the Stromgren's theorem, new solution of stellar structure equation can be obtained 
from the known ones through the homologous transformation. The new solutions describe 
new configurations with different masses, radius and chemical compositions (the so-called 
homologous stars). For recently obtained results see [s^ls^ where this problem is addressed 
in the context of brane cosmology. We pointed out that our Universe is also homological 
provided that it is filled by matter with the form of equation of state for noninteracting 
mixture of fluids p = —A + Cip(a) + C 2 a 8 and p = A + p ia~ 3( - 1+w ^ + p^a 5 + Po3 a ~ 3 - This 
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form is commonly used in cosmological considerations. The key idea of this work is deriva- 
tion this form of the equation of state from the first principles. It is proposed to derivate 
P = p{p) from postulate of self-similarity of its dynamics. This property is well known to 
engineers who build ship models as a prototype of real ship. It is justified by the fact that 
Navier-Stoke's equation are invariant against similarity symmetries. 

As a result we obtain the commonly used form of equation of state for mixture of nonin- 
teracting, dust matter, dark energy (the cosmological constant) and cosmic time variation of 
cosmological constant parameterised by the scale factor. Therefore subliminal role of math- 
ematics can also be seen in cosmology when we are looking for adequate form of equation 
of state for dark energy. 

We showed that the model with scaling multifluids fits well the supernovae data. For 
simplicity of presentation we demonstrate this for the case C\ = poi = 0. Physically it 
means that the Universe is filled generally by dust matter, cosmological constant fluid and 
additional scaling fluid which comes from Cardassian modification of the FRW equation (or 
equivalently it it is scaling fluid describing phantom fields for example). 

For the scaling multifluids model with A ^ 0, we found that Q m> o — 0.4 with n ~ — 1 
which correspond to (hyper) phantom model. If we assume fi mi o = 0.3 from independent 
extragalactic estimation, then Qq,o is small while value of n close to zero is favoured and 
model becomes close to ACDM. With assumption that A = the model is equivalent to 



42| on the new Riess sample. In 



Cardassian models. We confirmed our previous results 
particular the observations favour high density Universe with n negative. Again, if we 
assume fi m! o = 0.3 then n ~ which correspond to ACDM model. 

^From Fig. 1 can be seen the m-z relation for our model and ACDM one. We observed 
that distant SNIa should be brighter (in our model) than in the ACDM model. What 
is interesting that the Hubble diagram for the model under consideration intersects the 
corresponding ACDM diagram. In such a way the supernovae on intermediate distant are 
fainter then expected in ACDM model. This predictions could be tested with the future 
supernovae data. 

Our results demonstrate the existence of alternative model to ACDM model in explana- 
tion SNIa data. Therefore it should be interesting to compare both models from the point 
of view Akaike information criterion. Our result show that scaling multifluids model with 
vanishing A significantly better fits data then ACDM model. 
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FIG. 1: Riess Gold sample, M=15.935 Residuals (in mag) between the Einstein-de Sitter model 
and four cases: the Einstein-de Sitter itself (zero line), the ACDM flat model (upper curve), the 
best-fitted model (upper-middle curve), O m) o = 0.43, £lc,o = 0.28, Qa,o = 0.29, n = —3.33, and 
the best-fitted A = model (lower-middle curve) fi m ,o = 0.49, £lc,o = 0.51, n = —1.37. 

Moreover the cosmological model filled by scaling fluid makes a step toward solvin g th e 



coincidence problem of the present value of dark matter and dark energy components [24 1 . 

To make the ultimate decision which model describes our Universe it is necessary to 
obtain the precise value of fl m fl from independent observations because ACDM model and 
scaling multifluid models predict different density of the Universe. 
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TABLE I: Results of the statistical analysis of the model obtained both for the Knop et al. and 
Riess et al. samples from the best fit with minimum \ 2 (denoted with BF) and from the likelihood 
method (denoted with L). The same analysis was repeated with fixed O m; o. 



sample Q m ,o ^c,o n ^A,o M x 2 method 

K6 0.64 0.58 -3.17 -0.22 -3.61 53.4 BF 

0.38 0.64 -2.10 0.72 -3.53 L 

0.05 0.03 2.60 0.92 -3.53 55.1 BF 

0.05 0.03 0.43 0.92 -3.51 L 

0.30 0.13 -3.33 0.57 -3.55 55.1 BF 

0.30 0.00 0.00 0.70 -3.52 L 

K3 0.45 0.98 -0.47 -0.43 -3.49 60.3 BF 

0.27 0.01 0.00 0.75 -3.47 L 

0.05 0.05 1.97 0.90 -3.48 60.4 BF 

0.05 0.03 0.40 0.92 -3.46 L 

0.30 0.13 -2.23 0.57 -3.49 60.5 BF 

0.30 0.00 0.00 0.69 -3.46 L 

Silver 0.44 0.32 -3.27 0.24 15.895 226.7 BF 

0.46 0.40 -1.10 0.12 15.915 L 

0.05 0.07 1.77 0.88 15.935 229.4 BF 

0.05 0.07 1.50 0.88 15.945 L 

0.30 0.10 -3.33 0.60 15.915 230.9 BF 

0.30 0.09 0.00 0.61 15.945 L 

Gold 0.43 0.28 -3.33 0.29 15.905 172.1 BF 

0.44 0.34 -0.90 0.20 15.935 L 

0.05 0.10 1.57 0.85 15.945 174.0 BF 

0.05 0.09 1.30 0.86 15.945 L 

0.30 0.07 -3.33 0.63 15.925 175.2 BF 

0.30 0.00 0.00 0.70 15.945 L 
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TABLE II: Results of the statistical analysis of the model with fixed o = obtained both for 
the Knop et al. and Riess et al. samples from the best fit with minimum x 2 (denoted with BF) and 
from the likelihood method (denoted with L). The same analysis was repeated with fixed £l m< o- 



sample Q m fi Qc,o n ^A,o M. x 2 method 

K6 0.54 0.46 -3.33 0. -3.60 53.5 BF 

0.52 0.48 -3.33 0. -3.55 L 

0.05 0.95 0.33 0. -3.51 56.3 BF 

0.05 0.95 0.30 0. -3.51 L 

0.30 0.70 -0.10 0. -3.52 55.6 BF 

0.30 0.70 -0.13 0. -3.53 L 

K3 0.42 0.58 -0.77 0. -3.49 60.3 BF 

0.48 0.52 -0.40 0. -3.49 L 

0.05 0.95 0.30 0. -3.46 61.5 BF 

0.05 0.95 0.30 0. -3.46 L 

0.30 0.70 -0.13 0. -3.47 60.6 BF 

0.30 0.70 -0.17 0. -3.48 L 

Silver 0.50 0.50 -1.73 0. 15.905 227.1 BF 

0.51 0.49 -1.60 0. 15.905 L 

0.05 0.95 0.40 0. 15.975 239.3 BF 

0.05 0.95 0.40 0. 15.965 L 

0.30 0.70 -0.07 0. 15.945 232.3 BF 

0.30 0.70 -0.07 0. 15.945 L 

Gold 0.49 0.51 -1.37 0. 15.915 172.5 BF 

0.51 0.49 -1.23 0. 15.915 L 

0.05 0.95 0.43 0. 15.975 180.8 BF 

0.05 0.95 0.43 0. 15.975 L 

0.30 0.70 -0.03 0. 15.945 175.9 BF 

0.30 0.70 -0.03 0. 15.945 L 
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TABLE III: Model parameter values obtained from the minimisation procedure carried out with 
the Knop et al. and Riess et al. samples. 



sample ^ m ,o ^c,o ^A.o n 

K6 0.38i° : l 0.64±° : H 0.72lg;|° -2.10^;1 

0.05 0.031°;^ 0.92^;1 0.43l^ 7 
0.3 O.OOl ,;! 0.70+^;^ O.OOl?;^ 
K3 OJ27±g : S 0.0ll°;l 0.75±g^g 0.00±?;» 

0.05 0.921™ 0.40±S : gg 

0.3 0.00l°J8 0.69lg;» 0.00j3;» 
Silver 0.46±g:Sg 0.4O»% 0.12±gf2 -l-l^ 
0.05 0.07i°;°4 0.88l°£ l-SOlgj? 
0.3 0.09l°;£ 0.61±S;?? 0.001° |8 
Gold 0.44±g;il 0.34+o : | . 2 o+°jf -0.90l°;jj° 
0.05 0.09±g;^ 0.86l°;?l 1.30±g;g 
0.3 0.001°;?? 0.70lg;^ 0.001?$ 

TABLE IV: Model parameter values obtained from the minimisation procedure carried out on the 
Knop and Riess samples. Model with fixed Oa,o = 0. 



sample ^ m ,o ^c,o ^A,o ti 

+2.00 
-0.00 

+0.77 
-1.24 

+0.74 
-1.10 

+0.73 
-1.23 



K6 


n c 9 +0.09 
u - oz -0.09 


n 4 o+0.09 


0. 


-3.33 


K3 


0.4818;?* 


u ' oz -0.08 


0. 


-0.40 


Silver 




n 49+ - 05 

u '^ y -0.04 


0. 


-1.60 


Gold 




n 49+0.05 

u -^ y -0.04 


0. 


-1.23 



TABLE V: The Akaike information criterion (AIC) for models: ACDM model (ACDM), Cardassian 
and Scaling Multifluids. 



sample 


ACDM Cardassian 


Scaling Multifluids 


K6 


59.8 


59.5 


61.4 


K3 


64.3 


66.3 


68.3 


Silver 


236.6 


233.1 


234.7 


Gold 


179.9 


178.5 


180.1 
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FIG. 2: The density distribution for f2 m o- We obtain that m) o = 0.44^q'JJ on the confidence level 
68.3% (the inner dash lines). Additionally the confidence level 95.4% is marked (the outer dash 
lines). 
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FIG. 3: The density distribution for Q c ,o- We obtain that O, c ,o = 0.34I^§ on the confidence level 
68.3% (the inner dash lines). Both positive and negative values of Qco are formally possible. 



26 



TABLE VI: The Bayesian information criterion (BIC) for models: ACDM model (ACDM), Car- 
dassian and Scaling Multifluids. 

sample ACDM Cardassian Scaling Multinuids 
K6 63.8 65.5 69.4 
K3 68.3 72.3 76.3 
Silver 243.0 242.8 247.6 
Gold 186.0 186.7 196.3 
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„ 000 1 1 1 1 

-3.33 -2.33 -1.33 -0.33 0.67 1.67 
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FIG. 4: The density distribution for n. We obtain that n = — 0.90_ 1 90 on the confidence level 
68.3% (the inner dash lines). Additionally the confidence level 95.4% is marked (the outer dash 
lines). Both positive and negative values of n are formally possible. 
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FIG. 5: The density distribution for fi Aj0 . We obtain that Q Afi = 0.20 j^f on the confidence level 
68.3% (the inner dash lines). Additionally the confidence level 95.4% is marked (the outer dash 
lines) . 
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FIG. 6: Confidence levels on the plane (fi m ,o>^C,o) minimised over remaining model parameter. 
The figure shows of the preferred value of S7 mj o an d ^c,o- 



Tsujikawa S., Sami M., hep-th/0409212 



Gumjudpai B., Naskar T., Sami M., Tsujikawa S. hep-th/0502191 



Amendola L., Tocchini- Valentini D. astro-ph /001 1 243] 



Sahni V. Class. Quantum Grav. 19 3435, 2002 
Amendola L., |astro-ph7"9 904120 

Amendola L., Gasparini M., Piazza F. astro-ph/0407573 



Majerrotto E., Sapone D., Amendola L. |astro -ph/0410543 
Ellis G.F.R., Buchert T., |gr-qc/0506106| 



Harko T., Mak M. K. Phys. Rev. D69, 064020, 2004 

Copeland E.J., Liddle A.R., Wands D. Phys. Rev. D57, 4686, 1998 



29 




FIG. 7: Confidence levels on the plane (fi mi o, n) minimised over remaining model parameter. The 
figure shows of the preferred value of f2 m o and n. 



[32 
[33 
[34 

[35 
[36 
[37 
[38 
[39 
[40 
[41 
[42 
[43 



Falle S. Mon. Not. R. Astr. Soc. 250, 581 1991), 

Ibragimov N. K. Transformations in Mathematical Physics Moscow 1983 
Hydon P. E. Symmetry Methods for Differential Equation, Cambridge University Press, Cam- 
bridge, 1999 

Biesiada M., Szydlowski M., Szczesny J. Acta Cosmologica XVI, 115, 1989 

Lue A., Starkman G. D., Phys. Rev. Lett. 92, 131102, 2004 

Jafarizadeh M. A., et al. Phys. Rev. D60, 063514, 1999 

Silveira V., Waga I. Phys. Rev. D56, 4625, 1997 

Alam U., Sahni V., Starobinsky A. A. JCAP 04, 002, 2003 



Gorini V., Kamenschik A., Morchella U., Pasquier V., |gr-qc/0 403062 2004 
Godlowski W., Szydlowski M., Gen. Relat. Grav. 35, 2171, 2003 
Godlowski W., Szydlowski M., Krawiec A., Astrophys. J 605, 599, 2004 
Godlowski W., Szydlowski M., Gen. Relat. Grav. 36, 767, 2004 

30 



0. 040- 

c 
o 

D 

_q 0. 030- 
L 

-JJ 0. 020- 
>, 



0. 01 



c 

T3 



0. 000^^ 
- 1 . 00 




1 . 00 



C, 



FIG. 8: The density distribution for £lc,0 f° r the model with Q m) o=0.3. We obtain that £lc,o 



0.00 



+0.26 
-0.17 



on the confidence level 68.3% (the inner dash lines). Additionally the confidence level 



95.4% is marked (the outer dash lines). Both positive and negative values of £lc,o are formally 
possible. 



[44 
[45 
[46 
[47 
[48 
[49 
[50 
[51 
[52 
[53 
[54 
[55 



Dabrowski M.P., Godlowski W., Szydlowski M. Int. J. Mod. Phys. D13, 1669, 2004 
Godlowski W., Stelmach J., Szydlowski M. Class. Quantum Grav., 21, 3953, 2004 
Biesiada M., Godlowski W., Szydlowski M., Astrophys. J. 622 28, 2005 
Molina-Paris C, Visser M. Phys. Lett B 455 90 1999 
Tippett B.K., Lake K. |gr-qc/040"9088| 

Szydlowski M., Godlowski W., Krawiec A., Golbiak J. Phys Rev D72 063504, 2005 

Boyowald M. Phys. Rev Lett. 86, 5227, 2001 

Boyowald M. Phys. Rev Lett. 89, 261301, 2002 

Nojiri S., Odintsov S.D. Phys. Lett. B595, 1, 2004 

Elizalde E., Nojiri S., Odintsov S.D. Phys. Rev. D70, 043539, 2004 

Nojiri S., Odintsov S.D. Phys. Rev. D70, 103522, 2004 

Tolman R.C, Relativity Thermodynamics and Cosmology Oxford University Press, Oxford, 
1934 



31 




FIG. 9: The density distribution for n for the model with O mj o=0.3. We obtain that n = O.OO^'gg 
on the confidence level 68.3% (the inner dash lines). Additionally the confidence level 95.4% is 
marked (the outer dash lines). 

[56] Robertson H.P. Rev. Mod. Phys.5, 62 1933 

[57] Einstein A. Berl. Ber. 235 1931 

[58] Tolman R.C. Phys. Rev. 38 1758 1931 

[59] Steinhard P.J., Turok N. Nucl. Phys. Proc. Suppl. 124 38 2003 

[60] Shtanov Y, Sahni V. Class. Quantum Grav. 19, L101 2002 

[61] Coule D.H., Class. Quantum Grav. 22, R125 2005 

[62] Salim J. M., Perez Bergliaffa S. E., Souza ., Class. Quantum Grav. 22 975 2005 

[63] Pinto Neto ., Int. J. Mod. Phys. D13, 1419 2004 

[64] Barrow J. D., Kimberly D., Magueijo J., Class. Quantum Grav. 21 4289 2004 

[65] Bennett C. L. et al. Astrophys. J. Suppl. 148 1, 2003 

[66] Barris B. J. et al.Astrophys. J 602, 571, 2004. 

[67] Choudhury T.R., Padmanabhan T. Astron. Astrophys. 429, 807, 2005 

[68] Knop R. A. et al. Astrophys.J 598, 102, 2003 

[69] Riess A., Nugent P.E., Gilliland R.L., et al., Astrophys. J. 560, 49, 2001 



32 



0. ZAZ 

a 
O 

+-> 

D 

_q 0. 030 
L 
CO 

0. 020 

>, 



0. 01 



C 

T3 



0. 000^^ 
-0. 30 




1 . 70 



A, 



FIG. 10: The density distribution for Qa,o f° r the model with J7 mi o=0.3. We obtain that f2A,o 



0.70_ ' 26 on the confidence level 68.3% (the inner dash lines). Additionally the confidence level 
95.4% is marked (the outer dash lines). 



[70 
[71 
[72 
[73 
[74 
[75 
[76 
[77 
[78 
[79 
[80 
[81 



Riess A. G. et al. Astrophys. J. 607, 665, 2004 
Tonry J. L., et al., Astrophys. J. 594, 1, 2003 



Williams B. F. et al., 2003 astro-ph/0310432 



Freese K., Lewis, M., Phys. Lett. B 540, 1, 2002 
Zhu Z.-H., Fujimoto M.-K., Astrophys. J. 585, 52, 2003 
Sen S., Sen A. A. Astrophys. J. 588, 1, 2003 
Liddle A.R., |astro-ph/0401198| 2004 
Akaike H., IEEE Trans. Auto Control 19, 716, 1974 
Schwarz G., Annals of Statistics 5 461, 1978 

Parkinson D., Tsujikawa S., Basset B., Amendola L. Phys Rev D71 063524, 2005 
Jeffreys H., Theory of Probability, 3rd Edition, Oxford University Press, Oxford, 1961 
Mukherjee S., Feigelson E. D., Babu G. J., Murtagh F., Fraley C, Raftery A., Astrophys. 
J. 508 314, 1998 
[82] Stromgren B. Exapt. Ery. Naturw. 16 465, 1937 



33 



[83] Szydlowski M., Maciejewski A. J. Phys. bf A37, 1, 2004 



[84] http:/ / www-supernova.lbl.gov, http:/ / snfactory.lbl.gov 



34 



